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Abstract 

We examine the Markov properties of the three velocity components of a 
turbulent flow generated by a DNS simulation of the flow around an airfoil 
section. The spectral element code Nektar has been used to generate a well 
resolved flow field around an fx79w-151a airfoil profile at a Reynolds number 
of Re=5000 and an angle of attack of a = 12°. Due to a homogeneous ge- 
ometry in the spanwise direction, a Fourier expansion has been used for the 
third dimension of the simulation. 

In the wake of the profile the flow field shows a von Karman street like be- 
havior with the vortices decaying in the wake which trigger a turbulent field. 
Time series of the 3D flow field were extracted from the flow at different loca- 
tions to analyze the stochastic features. In particular the existence of Markov 
properties in the flow have been shown for different cases in the surrounding 
of the airfoil. This is of basic interest as it indicates that fine structures of 
turbulence can be replaced by stochastic processes. Turbulent and Marko- 
vian scales are being determined in the turbulent field and limits of standard 
Gaussian Langevin processes are being determined by the reconstruction of 
a flow field in time and space. 

1 Introduction 



Although computational methods for fluid dynamics have made enormous 
progress during the recent years, the problem remains to resolve fine struc- 
tures of highly turbulent flow. New turbulence models have improved the 
turbulence simulations to quite some extend [I][2]- Nevertheless, even with 
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the application of new techniques, it remains difficult to correctly simulate 
turbulent loads e.g. for wind turbines [3J[4J. Already the accurate simulation 
of the flow around single airfoils at higher Reynolds numbers is a difficult 
task [5] [6] as it is difficult to grasp the correct turbulence properties every- 
where in the flow. Therefore a lot of effort has been made to improve the 
modeling of the turbulence in small scales with subgrid models. The most 
common models are so called "deterministic" models based on eddy viscosi- 
ties. Meneveau et. al. give a good overview on such models for large eddy 
simulations (LES) [?]. 

Since a deterministic description of turbulence still has its limits, it is very 
convenient to describe the turbulence a time series [8] . The reconstruction of 
stochastic time series has made quite some progress in the recent years. First 
approaches involved the Fokker-Planck equation with Gaussian diffusion [9]. 
However today the research expands to Langevin equations for Levy pro- 
cesses [10] or even non-Markovian fields [UJ. From this evolved approaches 
to model turbulent data in different fields of science. Hence such approaches 
arouse some interest for the modeling of turbulent flows in numerical flow 
simulations \1'2\ . 

Some approaches to use stochastic models for flow simulations have been un- 
dertaken. Laval and Dubrulle proposed a Langevin approach for LES models 
which showed promising first results. Based on the Rapid Distortion The- 
ory (RDT)-model, described in [J3] [2] and [T5] , they developed a Langevin 
equation for an LES model, with friction made of viscosity and rapid distor- 
tion by resolved scales, using stochastic forcing with a mean value generated 
from the energy cascade. The model was validated on the flow in an empty 
cube against a DNS simulation. 

Mostly however, stochastic models were developed and used in the con- 
text of particel tracking or a Lagrangian framework as described by Bakosi, 
Pope, Shotorban or Fox [16] [T7J [18J [19] [20]. The particle tracking method by 
Bakosi works with a finite element grid in a Eulerian framework in contrast 
to the work of other groups. The physical magnitudes of the particles are 
described by probability density functions (pdfs) at each specific point. 
This approach requires the knowledge of the pdfs within a flow. One way to 
obtain these pdfs is the reconstruction by the Fokker-Planck equation using 
time series data. In an isotropic and homogeneous field of Gaussian turbu- 
lence this method is very easily applied and valid. In non-homogeneous flow 
fields the stochastic properties of the flow differ differ. A stochastic recon- 
struction of the fields using the method proposed by Friedrich and Peinke 
in [9] [ZJJ might nevertheless still be a promising approach to gain knowledge 
about the turbulent characteristics of a flow field for the modeling of pdfs. 
This shall be the focus of the following contribution. 
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2 The underlying equations for turbulence 
models 



The starting point for stochastic modeling is to consider a quantity in a 
statistical context. In fluid dynamics this is done using the Reynolds decom- 
position 

u = u + u' (1) 

where u is the averaged velocity field and u' denotes it's fluctuation. For 
time averaging the decomposition applied to the Navier-Stokes equations 
leads to the so called Reynolds Averaged Navier Stokes equations (RANS). 
For incompressible fluids they read: 

d t lli + djUjUi + dju'M = --dip + vd 2 Mi- (2) 

p J 

The incompressible continuity equation is then given as 

V • (u + vf) = V -u = V -u' = 0. (3) 
To solve the closure problem, models for the Reynolds stress term 

djWi ( 4 ) 

are needed. 

This is getting more complicated for LES models. Here the velocity field is 
spatially filtered and split into a resolved field Hi and a subgrid field given by 
u[ = Ui — Ui for Ui(xi,t) (note that here the average Ui describes an average 
over a spatial region and not - like in RANS - in time or over ensembles). In 
this situation the momentum equation of the Navier-Stokes equations turn 
out to be 



d t Ui + {uidi)uj + (uid^u'j + {u'fiijUj + dju'jiil = —d^p + vd]ui. (5) 

In this case a model for the forcing term / = (Hidi)u'j + (u^d^Uj is needed in 
addition to the subgrid scale stress tensor. Recently, Laval et al. [13] propose 
a stochastic approach to the first unknown term by putting up a Langevin 
equation: 

d t l =Al + £, (6) 
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where a Gaussian noise £ needs to be added in the equation, while A is a 
generalized evolution operator. Based on the observation that subfilter-scales 
are mostly dependent to resolved scales by a linear process similar to rapid 
distortion, Laval et al. derive from the momentum equation a model for the 
evolution operator A [T3] . 

Another approach is the one by Bakosi et.al. [IB] who model the stress terms 
using a (Gaussian) Wiener process to gain a pdf for the stresses. However the 
knowledge of the pdfs and the ability of their reconstruction are a premise 
for such a method. This requires the correct description of the turbulent field 
at all scales. Here the n-point correlation approach by using increments and 
the knowledge about Markov properties in a flow field is a basis [22] [23] from 
which further modeling can proceed. Tutkun et. al. have done an analysis 
of a spatial flow field experimentally |24j . 

As a typical example of a problem of turbulent flow we analyze here a flow 
field obtained by a DNS simulation of the flow around an airfoil to describe 
the method. The simulation, its parameters and the resulting field are pre- 
sented at first in the following section. Then the basic aspects of our stochas- 
tic approach is explained. Finally the results of the analysis of the data are 
discussed. 

3 The Simulation 

3.1 Simulation Parameters 

To gain a time series of a turbulent flow field a DNS simulation has been done 
using the high order spectral element code MenTar [25]. It combines finite 
element methods with the accuracy of spectral methods using Jacobi poly- 
nomials for the spectral expansion [26]. For the velocity pressure coupling 
a stiffly stable pressure correction scheme and an Adams-Bashfort predictor 
corrector time step were used as described in [27]. 

To obtain a well-studied and representative turbulent flow field from numer- 
ical flow simulations, the flow over an airfoil has been investigated [28]. The 
configuration is a 3D flow over a section of an fx79w-151a airfoil at an angle 
of attack of a = 12° at a Reynolds number of Re = 5000 in respect to the 
chord length. 

The used AfenTar code works with a 2D mesh with rather few elements for a 
first spatial discretisation. Further, the interior of the elements is calculated 
using an expansion by Jacobi polynomials [29] . Here we used a polynomial or- 
der of n=9. The discretisation in the homogeneous spanwise direction of the 
airfoil was realized using the "Fourier" version of the code which calculates 
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the third - spanwise - dimension using a Fourier expansion in combination 
with periodic boundary conditions on the sides of the domain (for validations 
see e.g. [SUES]). 

Measured in chord lengths L c , the domain size expands from — 10L C to 10L C 
in the cross-flow direction and from — 6L C to 20L C streamwise direction. The 
domain size in spanwise direction was set to n-L c , which showed good results 
in a simulation on the flow around a cylinder at Re=3900 [30J . The extension 
in the spanwise direction seems to be sufficient, as integral scales in the tur- 
bulent region in all directions are one order smaller than the domain size (see 
L in tabled]). 64 Fourier planes have been used for the expansion in spanwise 
direction together with a 2D hybrid mesh consisting of 2116 elements (see 
figHD- 

To resolve the boundary layer flow a quadrilateral grid layer was created 




Figure 1: 2D grid of the airfoil setup, consisting of 1179 quadrilaterals and 
937 triangles. A close-up on the right shows the boundary resolution and the 
resolution of the tail. 

around the airfoil. The original shape of the airfoil was not changed since 
every detail of the geometry has been taken into account. As the tail of the 
airfoil had a tiny triangular shape, the resolution of this tail leads to very 
small elements at the trailing edge (see fig{T]on the right). This resulted in 
very small time steps of order t w 10~ 5 . In our simulation the maximum 
CFL was < 0.4. 

To avoid instabilities which in turbulent flows often develop from Neumann 
outflow boundary conditions a sponge layer has been implemented increasing 
the viscosity in the last 3L C before the outflow. Finally we used a uniform 
laminar velocity distribution at the inflow. 
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3.2 Simulation Results 

Fig. [2] depicts the average flow field over the airfoil. Here and further on, 
all units are dimensionless and normed to chord length of the airfoil. The 
average flow field is mainly characterized by a laminar flow separation in the 
tip region with an impingement point in the area behind |L C . The shear 
flow which develops downstream the separation point leads to a von Karman 
like vortex street in the rear part of the airfoil. Triggered by the vortices a 
turbulent wake flow field appeared (see e.g. vorticity contours in fig. EK))- 
To give an impression of the distribution of the turbulence intensity in the 
field, fig. [3]d) shows contours of the RMS of the velocity magnitude of the 
flow in the vicinity of the airfoil. In this figure high fluctuations in the veloc- 
ity magnitude are recognizable mainly close to the trailing edge and in the 
wake. Thus for the study of turbulence at and around the airfoil the regions 
at the trailing edge and in the wake of the airfoil were of main interest. 
To perform a time series analysis of the flow properties specific points effected 
by the turbulent wake flow have been chosen within selected elements of the 
mesh in the region depicted in fig. HI 

The regions selected were positioned on one hand in the wake (element 




Figure 2: Mean velocity of the flow field over t=20 in normalized time scales. 
The magnitudes of the velocity given in a) by the colors, additionally the 
vector orientation of the flow in b). 

1), on the other hand close to the trailing edge of the airfoil at different 
heights (elements 2-5). It must be noted that all position selected for the 
recording of time series are located in regions with a triangular unstructured 
grid. Within the selected elements at the spanwise position of z = 2.11L C a 
time series of 10922 datasets at 208 different spatial points were collected at 
a sampling rate of f r = 400. Out of the five selected elements 62 points have 
been chosen for the purpose of the following analysis. 
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Figure 3: Snapshot of vorticity components in z-direction a) and rms- 
fructuations of the velocity magnitude in b). 




Figure 4: Selected grid in the region just above the trailing edge for the 
further investigation of the turbulence 
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4 The Definition of Markov Properties 



In the following the stochastic analysis of the turbulent flow field will be 
discussed. The final aim is to gain more knowledge of the Reynolds stress 
tensor, which for a the stochastic analysis reads 

J wJ«j-p(«X-) d ( u i u j)> ( 7 ) 



u i u j 



where piu'^u'^) is the probability density function (pdf) of the components of 
the Reynolds stress tensor u'jUj . 

In the following we will generalize the discussion on the stochastic features 
of important fluctuating quantities such as Reynolds stress. For turbulence 
it is of central interest to characterize the spatial correlations of quantities 
4>(xi,t). Therefore it is common to investigate the statistics of increments in 
a scale of r at the same point in time t 

<j>(r,t) = <f)(xi-r,t) - <f>(xi,t), (8) 

where Xi denotes the position in space and t is the time. 

The following analysis will be based on probabilities of such increments at 

different scales. The conditioned probability function of can be described 

as 

/ \ p(0i,n;02,r 2 ; ...;<f> n ,r n ) 
p(<f>i,n\(f>2,r 2 ;...;<f) n ,r n ) = . , x (9) 

P{(P2, T 2, <Pn, T n ) 

where the left hand side reflects the probability density function of the mag- 
nitude 0i at the incremental scale r\ under the condition that at the scales 
r 2 ,...,r n the magnitudes 2 ,...0 n have selected fixed values. Here we use 
Ti < r i+ i and the abbreviation <f>i = <f>(ri,t). 

In stochastics we speak of Markov properties, if the pdf of a dataset can be 
completely described with the knowledge of the neighboring magnitude at 
scale r n _i without being influenced by the next magnitude at r„_ 2 - In case 
of the increments this means: 

p(0n, r„|0„_i, r n _i; fa, r\) = p(0„, r„|0„_i, r n _i) (10) 

If this is the case the general n-scale statistics can be closed by the two 
scale statistics p((f) n , r„|0„_i, r n _i). Furthermore, it is possible to apply the 
Kramers-Moyal expansion of the data series, which leads to: 

= ±(-^) k D«HMp(M M . (ii) 
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The coefficients of this expansion are defined as 

D {k) ((f),r) = lim M (fc) (0,r, Ar) (12) 



Ar^O 



This definition of D^n) can be used to estimate its value directly from the 
data [9] [21] [22]. 

The theorem of Pawula states that = V k > 2 if the coefficient of 
the expansion = [32]. In such the Kramers Moyal expansion 

ffTTj) truncates to the second order, yielding the Fokker-Planck equation (also 
called Kolmogorov forward equation) 



d rP (4>, r |0 O , r ) = [-fy£>Wp(0, r) + d^D^p(<P, r)) P (<P,r\(f> , r ). (14) 

The first coefficient .D^ is referred to as the drift coefficient, the second 
coefficient D^ 2 ^ is called diffusion coefficient and determines the strength of 
a Gaussian distributed white noise. This noise is 5-correlated due to the 
Markov properties [32J. The Fokker-Planck equation can be reformulated 
as a Langevin equation following Ito or Stratonovich (see [32J ) . The Ito 
formulation is: 



°" {r] - 1 -D (1) (0,r) + J-D( 2 )(<P,r)T(r). (15) 



v f 

T(r) is the Gaussian distributed white noise. 

Under these conditions it is possible to analyze the Reynolds stress tensor 
4> = u'jUj. Having determined and D^ 2 \ the Fokker-Planck eq. (fT4"j) 
allows to reconstruct a series on all n-scales of r n [33]. Thus it is possible 
to reconstruct p(0j,rj) for given scales r^, where may be given by a finite 
spatial relation. This can be used as an expression for the components of 
the Reynolds stress tensor u^u'^r) eq. (jlj). A similar procedure could be 
constructed for other relevant turbulent quantities like the additional stress 
terms in equation ©. 

The aim is now to find out under which conditions this stochastic simplifi- 
cation holds for a given flow problem. However, the geometrical distribution 
of the data points of the flow field makes a straight forward analysis difficult 
since not enough straight r n scales could be evaluated for a complete recon- 
struction. Therefore the analysis is being done using the Taylor hypothesis 
of frozen turbulence, which is applied quite frequently also in experimental 
works (like e.g. [22] [23]), by 

r = tu. (16) 
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Now, r = (t— to) is a time scale by which the same analysis can be undertaken 
as in given above in equations flU])-f lT5]) . In this paper we focus on the most 
important underlying Markov property Thus we look if 

p(u , T \u h n; u 2 , r 2 )) = p{u , T a \u U Ti) (17) 

in time and in spatial scales if 

p(u , r \ui, ri; u 2 , r 2 ) = p(tt , r |«i, r x ) (18) 

holds, u is one velocity component and r = xq — x\ is a spatial scale for 
this coordinate. For simplicity reasons the nomenklatura used for the coor- 
dinates and orientations has been chosen as in usual Cartesian coordinates 
giving (x, y, z) for the streamwise, crossflow and spanwise direction as well 
as (u(t),v(t),w(t)) for the velocity increments in this direction. 



5 Analysis of the results 

From the over all flow field four points (marked as 13,61,119 and 218) have 
been selected for further evaluation (see fig E]). In table [1] the scale of the 
spacing between the points r is given for all components (x,y,z) as well as 
the mean velocity components ul, standard deviation a, of the velocities and 
estimations of the integral, Taylor and dissipation length scales, represented 
by Li, Xi and rji respectively. All length scales have been calculated accord- 
ing to [3l]. Since the Taylor length is a fit with possible deviations and the 
dissipation scale is derived from the Taylor length, both more or less give 
values for the order of magnitude and are not to be seen as exact values. 



5.1 Analysis in time scales 

The first analysis of the time series has been conducted using the dataset of 
the point marked with the number 218 as a reference in time scales r. We are 
later interested to use the Taylor hypothesis. The mean flow field direction 
is oriented in the x-y-plane. Therefore we focussed here on the analysis of 
the u- and v-velocity increments. However, the complete analysis has also 
been done for the spanwise direction, even though for brevity these results 
might not all be presented here. 
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Figure 5: Geometric positioning of the evaluated points in the region over 
the tail of the airfoil. The numbered points have been selected for further 
analysis and reconstruction 

To give a first idea of the distributions of the velocity increments, the non con- 
ditioned histograms of the single components are shown in figJS]for r = 0.15. 
It is obvious, that for the u and v-components there is some skewness in the 
distribution and all distributions are non-Gaussian. 



Next we discuss at which scales the flow field does show Markov properties. 
For all geometrical points a dataset of a time series of 10922 points in time 
was recorded. Since this was very short, only 15 bins have been used for the 
evaluation. 

FigJT] shows the contours of the conditioned probability density function for 
p(ui,ToK,i,Ti) and p(ui,T \ui A ,T 1 ;Ui^,r 2 ) for At = t 1 -t = 0.15, r = 0.15 
and T2 = ro + 2Ar for the velocity components u and v. Also slices of the 
contour-plot are given to give an impression of the shape of the conditioned 
pdfs. The shapes of the pdfs do not appear very smooth as it has to be 
kept in mind, that the over all number of data is just above 10000 and the 
plots are logarithmic. However, the similarity of the two contours suggest 
the existence of Markov properties. 



Next, a more quantifying test is done to check if Markov properties are given. 
To evaluate the Markov properties, now a x 2 -test has been performed for the 
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Figure 6: Histogram of the velocity increments in semi- logarithmic presen- 
tation for u in a), v in b) and w in c) and r = 0.15. 
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Figure 7: a) and c) give the contour plots of the pdfs of p(ui, ro|«i,i, Ti) 
(gray) and p(ui, ro|wi,i, T\\ u it2 , r 2 ) (black) for r = 0.15 for u- and v-increments 
respectively. Slices at about +/— 0.8a for u- and +/— 0.9a for v-increments in 
b) and d) respectively show a high coherence between the two distributions, 
where the dotted line represents the data from p(ui, TqIm^i, Ti, u^, T2) and 
the errorbars are given by y~N with the number of events in the bin. 
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Table 1: The analyzed points, their geometrical distances - with point 218 as 
a reference point -, their velocity and the statistical properties of the velocity 
for all components. 



distributions of p(ui, Tol^i, 7~i) and P2{u i) r\ui$, r ; m^i, Ti) (see [35J) 

r = E — 7f f - («) 

Here pi(«) represents the first probability density - in this first casep(wj, r \ui t i, Ti) 
i is the bin number and P2{i) is the second probability density respectively 
- here p(ui, Tolw^i, t\\ Ui t 2, t 2 ). The test quantifies the difference between two 
pdfs from data with errors. Usually as p\ a theoretically know pdf is taken. 
For our purpose of the estimation of the existence of Markov properties, this 
might lead to cases where p\{%) = 0. In such cases we set x 2 has been set to 
0. Thus x 2 — events had to be evaluated separately. 

x 15 bins have been used for evaluation, giving 15 degrees of freedom. Taking 
the average of the x 2 values for all the evaluated distributions shown in fig|8] 
for all velocity components, we obtain an impression of the validity of the 
assumption of Markov properties. 

In figJH] it becomes evident, that the bins containing few incidences fur- 
ther from the mean values, show stronger deviances in the distributions. 
This might be due to statistical reasons which lead to very low binning val- 
ues, however it might also be caused by a non-Markovian behavior of the 
stronger fluctuations in the flow at some scales. This could only be clarified 
if the dataset would be increased for a longer simulation period. 
Due to this insecurity in the Markov analysis, five bins have been excluded 
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Figure 8: x 2 values for the conditioned pdfs with r and Ar = 0.15 against 
the u\, Vi, and w\- velocity increments bin. x 2 is calculated as the sum over 
all velocity bins u, v and w. The line marks the value under which the 
probability of the distribution to be non-Markovian is less then 5% in the 
classical \ 2 theory. 



from the evaluation for the further analysis: The three outer bins to the 
lowest negative velocity increments and the two bins containing the highest 
velocity increments. FigJH] shows the dependency of the average x 2 on Ar 
for 0.0125 < Ar < 0.45 for the velocity components, a) and c) including all 
bins and in b) and d) using only 10 bins after the exclusion of the five bins 
containing to few data. Note, that with the decrease of degrees of freedom, 
the limit for the criterion of the probability being < 5% for non-Markov 
properties also decreases. 



Using the 10 bins we can expect Markov properties at Ar > 0.15 for the v- 
and w-component while for the u-component we supposed to have a Marko- 
vian field for 0.075 < Ar < 0.375. Applying the Taylor hypothesis (eq. 
(TT5|) ). this would correspond to a so called Markov length of Xm = 0.011 
for the u-velocity and Am = 0.018 for the v-velocity respectively. Note the 
comparable values to A in table [1] as proposed by [36J . 

Within the estimated validity range for the Markow properties, now the 
Kramers-Moyal coefficients have been determined for each r and Ar. There- 
fore for each r a fit for ZimA T - >o corresponding to equation fflBl has been 
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done. To estimate a linear fit was used, whereas for and log- 
arithmic fits have been used, since they fitted to course of the curve better 
and avoided unphysical negative values. As an example in fig HQ] a), b) and 
c) the fit for the normed v- velocity increments at r = 0.15 is shown for 
and DW. 



The D^(ui,r) and D^ 2 '(ui,r) have further been fitted by a linear and a 
quadratic function over all described by Renner in [22J for isotropic 

turbulence (see fig {10] d) and e)). These fits were in the following used for 
the reconstruction. 

Even though we cannot assume D^(u{,t) = in the dataset used, since 
« we take the Fokker-Planck equation (JT3J) as an approximation. 



5.2 Reconstruction in time 

To get an estimation of the accuracy of a reconstruction assuming a Gaus- 
sian field, the reconstruction has been done using only the first two Kramers- 
Moyal coefficients. 

The obtained D^(ui,r) and D^ 2 \u^t) are used to reconstruct the time se- 
ries with a Langevin equation ()15p modified for time scales. In figOJ] the 
pdfs p(ui, ro|wi,i, Ti) and p(ui, ro|u»,i, T\\ Ui irec , r 2 ) of the u- and v-velocity in- 
crements for t = 0.15 and T\ = 0.30 are compared. Here u itrec (T 2 ) is the 
reconstructed data at a 3r scale. Below again slices at about +/ — 0.8a for 
u and +/ — 0.9cr are shown. 



The pdfs of the increments show a non perfect, but quite good reconstruc- 
tion. Fig. [12] shows the semi-logarithmic plot of the pdf of the original data 
against the reconstructed data for r = 0.15. Although certainly the dataset 
was still short, the general shape of the curve is met in most cases (see table 
15.21) . The main reason for deviations is most likely the approximation for the 
Kramers-Moyal expansion, since ^ as the higher order Moments are 
not met perfectly. 
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Table 2: Statistical properties of the distributions of the original data 
(marked with data) and the reconstructed time series (recon) for r = 0.15. 
Notable is that the higher order moments are not reconstructed correctly. 



The contours of the single conditioned pdfs p(u, Tq\ui, ri) and p(u, T \u rec , Ti) 
give a good impression of the quality of the reconstruction (figHSJ) in the 
relation between the developments of the pdfs in time. While the overall 
tendency of the structure of the main peaks and their form is very similar, 
the spatial extension differs a bit between the two pdfs. 

5.3 Reconstructing the spatial distributions 

One of the main problems for flow simulations is however the unknown flow 
statistics of the flow at a distant point. After having obtained convincing 
results for the Markov properties of velocity increments, we will apply next 
the Taylor hypothesis to investigate the velocity increments in for the selected 
points 218,61,13, and 119 of fig. El 

For local isotropic turbulence it is well know that for increment statistics in 
space and in time are related by the Taylor hypothesis. Here we will consider 
a flow with a shear which is not homogeneous or local isotropic. The open 
question is, whether it is possible to use the above obtained characteristics 
of the turbulent flow for an application in space. 

First we start with the Markov properties. Fig. O shows the contours of 
the conditioned pdfs p(«j, tqIm^i, r x ) and p(ui, r \u it i, r±; u^, r 2 ) of the u- and 
v-velocity increments between the points 218-61, 218-13 and 218-119, with 
overall distances r = 0.054, ri = 0.076 and r 2 = 0.088. 



The contours in figHH indicate good Markov properties. This is reflected in 
the x 2 values for the two distributions shown in figE] where some of the 
outer bins have a lack of data points. However, especially close to the mean 



17 



velocity increments Ui ro the Markov properties can be very well assumed as 
the probability of a non Markovian field is < 5%. 

Taking the result as a strong indication that Markov properties hold we 
want to show that the spatial statistics can be reproduced by the stochastic 
processes we extracted in the previous section from the increments using a 
time step r. We reconstruct from the results obtained for the point 218 data 
series using eq. (1151) . For a r = - where r = d(218, 119) the correspondig 
conditional pdf can be obtained for all components (see fig JT6|) . 



Further we compare the pdfs of the increments of r 2 = d(218, 119) of the 
simulated data with the reconstructed increments. The results for all com- 
ponents are given in figJT71 As can been seen the spanwise flow direction 
was obviously the easiest to reconstruct. Since there was no big difference 
velocity magnitude between the two points and the original distribution was 
also symmetrical. For the pdfs of the u- and v- velocity increments show good 
agreement, even though slight deviations can be observed. 



By adding the reconstructed increment to the velocity time series, one can 
finally resolve a reconstructed time series for the velocity at point 119: 

with Uiiy&if) being the Mj- velocity component at the point 218 at the time 
t and u r 2,rec the reconstructed increment for r 2 = d(218, 119). Thus we are 
able to compare the characteristics of the original velocity at point 119 given 
in table [I] to a velocity reconstruction. Table 15.31 displays the results of this 
reconstruction. Concerning the mean, the u-velocity shows some deviation. 
Although the tendency to reduce the velocity closer to the airfoil surface 
is correct, this reduction is not yet large enough. On the other hand the 
values for the v- and w-velocity are met quite well. This contrast leads 
to the conclusion, that mainly the u-component had effects in the flow field 
which were not grasp by the Gaussian statistics. Nevertheless the results was 
surprisingly good as in an inhomogeneous, anisotropic field with a physical 
object in the area it was possible to reconstruct the field by pure statistical 
analysis. 
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Point 


Component 


v~ 




U 


\ 


Vi 


119-rec 


X 


-0.093 


0.292 


0.162 


0.0015 


0.0001 


119-rec 


y 


-0.042 


0.199 


0.127 


0.0005 


9.4e-05 


119-rec 


z 


-0.002 


0.192 


0.115 


0.0002 


6.8e-05 



Table 3: Length scales, mean and standard deviation of the reconstructed 
time series at point 119. 

6 Conclusions 

We analyzed a simulated flow field with the aim to gain information about 
the statistical properties of the field for the reconstruction of the statistics 
of the flow as it might be useful for stochastic models in computational fluid 
dynamics. Further the attempt was made to reconstruct the statistical prop- 
erty of the flow at one point from the stochastics at a different point. The 
field was obtained by a simulation around an airfoil using the spectral ele- 
ment code MekTolt. 

The simulation generated a well resolved flow field with an anisotropic and 
inhomogeneous turbulence area in the wake of the airfoil. The resolution 
of the field enabled us to use the flow field for an analysis of the stochastic 
characteristics of the turbulent wake using n-point correlation methods with 
velocity increments. 

In a first step the time series has been analyzed at a point within an anisotropic 
and inhomogeneous turbulent field in the wake of the airfoil. Although the 
number of data points was low there was a strong indication that there exist 
Markov properties for time scales greater than a certain "Markov timescale" . 
Applying the Taylor hypothesis this would correspond to a Markov length of 
Xm = 0.011 for the u-velocity component and Am = 0.018 for the v-velocity 
component respectively. However due to the lack of data an uncertainty on 
Markov properties remains for the more rare incidences of large increments. 
Assuming that Markov properties hold, the conditional moments M^ n \u, r, At) 
and Kramers-Moyal coefficients D^ n \u, r) for the time series for different At 
and r have been calculated. = was not given, however « 
was given. Hence the further analysis has been conducted under the as- 
sumption that the conditions of Pawulas theorem were met. The first two 
Kramers-Moyal coefficients were used for a reconstruction of the incremental 
time series at the analysis position. 

This reconstruction of the incremental time series by a Langevin equation 
gave promising results. Slight deviations were observed since we used Gaus- 
sian methods for reconstruction. Thus skewness and kurtosis did show some 
differences. 
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Further we used the obtained coefficients to estimate the structure of the flow 
field in real space of separated points within the inhomogeneous turbulent 
field by applying the Taylor hypothesis. Therefore the incremental analysis 
was also done using four points aligned in the wake including the original 
analysis point. The distance between the base point and the position of the 
reconstruction was about 5-8 times the Markov length scale Xm- 
The results were surprisingly good, especially for the v- velocity component. 
As the flow field is very inhomogeneous due to the shear flow on the one 
side of the field and a boundary layer towards the airfoil, the reconstruction 
just based on the statistics of the field in one point using a basic Langevin 
equation, the results quite good. 

Since the results were quite promising, we do encourage further research 
in this area. It would be especially promising to extend the model to non 
Gaussian Langevin equations for the inclusion of the effect of higher order 
moments. It also should be investigated to which accuracy and on which 
scales it is possible to assume the statistics of a distant point from the statis- 
tics in another point - within an inhomogeneous filed. The results in this 
work give hope that such methods could lead to quite satisfying results, par- 
ticularly since the dataset was still very short for the used method. 
Towards the aim of implementing the method into a general stochastic turbu- 
lence model also further questions need to be answered, such as how the flow 
properties change in the boundary layer, the direct shear area or for different 
Reynolds numbers. Answering these questions should make the method im- 
plementable in stochastic models for flow simulations as the ones proposed 
by Laval or Bakosi [13j[16j and increase their accuracy. 
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Figure 9: x 2 in dependency of Ar for all bins a) and c) on the left and 
for only 10 well populated bins on the right b) and d), as a) and b) depict 
the u-component and c) and d) the v-component. It becomes evident, that 
Markov properties can be assumed in case of the smaller absolute value of 
the increments where enough data is available. The assumed validity of the 
Markov properties in the date here is Ar > 0.15 for the v-component and 
0.075 < Ar < 0.375 for the u-component. 
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Figure 10: Examples of the fits of M (1) , M (2) and M (4) in fig. a), b) and c) 
respectively for r = 0.15. The used range is plotted in gray, while black dots 
have been excluded. The resulting normed and plotted on all 

used bins with fits for D^> and D^> in d), e) and f). 



25 




-0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 
v1 v V 

Figure 11: Comparison of the contour plots in a) and c) of the pdfs 
ofp(ui,r \ui t i,Ti) and p(ui, ro|«i,i, Ti; Ui 7rec , T2) for r = 0.15 in case of the 
reconstructed time series Ui >rec for u- and v-increments. Again slices at 
about +/ — O.80" for u and +/ — 0.9cr for v in b) and d) respectively 
give an impression of the distributions. Here the dotted line represents 
the data from p(ui,To\ui t i,Ti] Ui^T^) and the solid line the data from 

p{u i ,T \u iA ,T 1 ]U^ rec ,T 2 ). 
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Figure 12: Logarithmic plot of the non conditioned pdfs p(u, r ) (line) and 
p(u rec , To) dots at r = T\ = 0.15 showing good general coincidence and some 
differences in particular due to the different skewness of the distributions. 




Figure 13: Contour plot of the single conditioned pdfs p(u,r\u ,T ) (black) 
and p(u,T\u rec ,To) (grey) at r = t = 0.15 showing the coincidences and 
some differences between the contours, especially due to the asymmetry of 
p(u,t\u ,t ). 
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Figure 14: Contour-plot of conditioned pdfs of the time series p(ui, r | uj,i, r%) 
(grey) and p(ui, ro|ui,i, r\\ «i,2, ^2) (black) for the u- and v-velocity increments 
between point r0=d(218,61), rl=d(218,13) and r2=d(218,119). 
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Figure 15: x 2 values for the conditioned pdfs for the increments between 
points 218, 61, 13 and 119 with overall scales r = 0.054, r = 0.076 and 
n = 0.088. 
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Figure 16: Contour-plot of conditioned pdfs of the time series p(ui, r |uj,i, ri) 
(grey) and p(ui, ro|tti,i, r±; Ui jrec , r-z) (black) for the v-increments a). Slices of 
the plot at +/ — 0.9cr for the v component. Here the dotted line repre- 
sents the data from p(ui, ro|«j,i, r±; Ui t 2, r-z) and the solid line the data from 
p(ui,r \ui ! i,ri;u ijrec ,r 2 ). The points with error bars are obtained from the 
data of the simulation with the errorbars being yN. 
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Figure 17: The semi-logarithmic plot of the pdf of the velocity increments for 
all components, where the line is the distribution coming from the original 
simulation data and the dots represent the reconstructed time series. 
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